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We present a number of second order maps, which pass 
the singularity confinement test commonly used to identify 
integrable discrete systems, but which nevertheless are non- 
integrable. As a more sensitive integrability test, we propose 
the analysis of the complexity ("algebraic entropy") of the 
map using the growth of the degree of its iterates: integra- 
bility is associated with polynomial growth while the generic 
growth is exponential for chaotic systems. 

05.50.+q,02.90.+p 



Discrete systems have for a long time been a subject of 
study in the field of dynamical systems. A strong prac- 
tical motivation is the power of numerical exploration 
for such systems Q, which led to interesting findings in 
chaos theory. On the other hand, several numerical al- 
gorithms ( "convergence acceleration algorithms" ) are re- 
lated to integrable discrete maps (see 0] and references 
therein). Finding out whether a given system is chaotic 
or integrable is then a basic question worth more investi- 
gation. One open problem to find an algorithmic test of 
integrability, and this is the subject of the present letter. 

In the case of continuous systems one can test for the 
"Painleve property" [[|-|| , which is closely related to in- 
tegrability and has considerable predictive power. An 
analogue of the Painleve test for discrete systems was 
proposed in j?]] and has been used as a powerful construc- 
tive tool, e.g., to identify discrete Painleve equations ||. 
This singularity confinement test is similar in spirit to 
the continuous Painleve test in that it analyses behavior 
around a movable singularity of the map. When a map is 
iterated, it may happen that we reach a point for which 
the next value is ill-defined due to the appearance of an 
indeterminate form, oo — oo, • oo or such. One should 
then study the behavior around the singularity: if the 
map can be continued in a way which allows, after a fi- 
nite number of steps, to exit from the singularity without 
loss of information, then the system is said to pass the 
test. 

We show in this letter that the confinement test is not 
sufficient to ensure integrability. We also propose another 
indicator for rational maps: a measure of the "algebraic 
entropy" |] , which has to do with global properties of 
the system (see later). 

We shall consider the map 

x n+ i + x n -x = x n + a/x 2 nl (1) 
and some of its generalizations. Relation dl|) defines a 



map (i„_i, x n ) — » x n+ \. The potential singularity of this 
map is reached if, at some step (say step 0) , we arrive at 
xq = (with a finite nonzero x—i), because then x\ = oo, 
X2 = oo, but X3 = 00 — 00 and it is not clear how to 
proceed. To refine the analysis, let us assume that we 
arrive at xo = e by x-\ — u, with suitable previous x^s. 
With these initial values we get the sequence (here a=l) 

X-l= u, 
x = e, 

x\ = e~ 2 — u + e, 



X-2 



e- 2 -u 



x 3 = -e + 2e 4 + . . . , 
X4 = u — e + 

In this case the outcome is: "00 — 00 = 0" and the se- 
quence emerges from the singularity with the value it, 
i.e., without losing the initial information. This means 
that the system (|l|) passes the singularity confinement 
test without problems; the singularity structure u — > e — > 
e~ 2 — > e~ 2 — > — e — * u is rather typical. 

The problem is that system ([!]) is chaotic, as we shall 
show. Our suspicion about the non-integrability of (Q) 
arose when, motivated by [[Tof , we evaluated the growth 
of the degree of its iterates as follows: 

We start by writing the map as a first order two- 
dimensional map 



V ■ Pn — (x„_i,x n ) — ► Pn+l = (x n ,x n+ i), 



(2) 



and then rewrite tp in terms of homogeneous coordinates 
[y n ,z n ,tn] by setting 
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(3) 



This means that we are now working in the two- 
dimensional projective space CP2, and that points with 
homogeneous coordinates [y, z, t] and [Ay, A z, At] are to 
be identified (projectivization). For ([!]) the map <p may 
be written as 



(4) 



In CP2 the above singularity pattern looks as follows: 





y 




y 3 + at 3 — y 2 z 


<p ■ 


z 




y 3 




t 




ty 2 








1 




1 












u 









1 




1 







1 























1 



The last term of this sequence is not in CP2 and is now the 
manifestation of the ambiguity mentioned above. Note 
also that in this formulation it is clear that infinities are 
not singularities: they look like any other point, the last 
component is just zero. The expansion around the sin- 
gularity clarifies the situation. We get the sequence 
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3ue 2 
Sue 2 



9ue 2 



u - 
-e ■ 
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and in the last term we are able to cancel the factor 
e 8 , after which we can let e — > 0, getting [u, 0, 1]. We 
have thus emerged from the singularity with the initial 
information u. 

The cancellation mentioned here is crucial. It occurs 
only if there is a singularity in the map, because the 
existence of a singularity means that there will be com- 
mon factors of e. Such cancellations of common factors 
are necessary to reduce the growth of the degree, because 
otherwise the successive iterates (p^ of tp would be poly- 
nomials of degree d n , where d is the degree of p. For 
integrable systems the cancellations are in fact so strong 
that asymptotically the degree grows only polynomially. 

The degree of <p> is not canonical, since it is not in- 
variant under coordinate changes. However the growth 
of the degree is canonical It is generically expo- 

nential fll|l , but may become polynomial if the number 
of common factors is large enough. The conjecture is 
that intcgrability of the map implies polynomial growth 
@(scc also fl|). 

For the map (||) we get the following sequence of de- 
grees: 1, 3, 9, 27, 73, 195, 513, 1347, 3529.... The 
first four degrees follow the 3™ rule; cancellations then 
take place. Note that the first drop of the degree is 
3 x 27 — 73 = 8 corresponding to the factorization of 
e 8 in the above calculation. From the nine first numbers 
of this sequence we inferred the generating function for 
the degrees: 



1 + 3^ 



(1 - x) (1 +x) (x 2 - 3x + 1)' 



(5) 



The next degree, found by iteration of ip is 9243 and 
coincides with the prediction obtained by expanding g: 

g = 1 + 3 x + 9 x 2 + 27 x 3 + 73 x 4 + 195 x 5 + 513 x 6 + 
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+1347 x ' + 3529 + 9243 x a + 24201 x 

(A proof that (bh indeed is the generating function of the 
degrees will be given elsewhere [pj.) 



Function (||) generates a sequence with exponential 
growth. The denominator of g contains two basic pieces 
of information. It first shows that the sequence of degrees 
verifies the very specific relation 



d ■ d. 



n+3 



d ■ d 



n+l 



0. 



(6) 



where d = d\ = 3 is the degree of the map ip, and d n 
the actual degree of ip n , after "projectivization" . It also 
determines the asymptotic behavior of d n : if a is the 
smallest modulus of the roots of the denominator of , 
then d n +i rs a~ 1 d n asymptotically. In this case a = 
(3 — \/5)/2 and we define the "algebraic entropy" of the 
map by 



£ = lim — log(c?„ 



log 



(7) 



This calculation indicates that the map has non- 
vanishing entropy, and therefore is likely to be non- 
integrable. 
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Fig. 1: A collection of orbits of the map (Qj 

Numerical chaos can be seen when we draw a picture 
of some orbits of the map, see Figure 1. This figure was 
obtained with a = 7. The two "cat's paws" are around 
two points of a nine-periodic orbit of ip. Such a point 
of order nine of <p (x n — X n +\ ~ 3.043896 . . .) is located 
at the lower left hand corner of Figure 2. The picture 
is characteristic of chaotic behavior of a two dimensional 
conservative system. 

In order to compare with a truly integrable system let 
us consider (T3l 



Xn+l ~T X n —\ 



a. 

~I2 



(8) 



which is related to d—Pj. This is also a third degree map, 
but the singularity structure is now such that the degrees 
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grow only as 1, 3, 9, 19, 33, 51, 73, 99, 129, 163, 201, ... The 
generating functional is 



and in fact the degrees grow polynomially according to 
the simple rule 

d n = 2n 2 + 1. 
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Fig. 2: A close up of Figure 1 

If model (Q) were an isolated example, it might be dis- 
missed by some ad hoc rule. However, it turns out that 
the singularity confinement test is somewhat insensitive 
in general, and maps containing rather arbitrary func- 
tions pass the test. Let us consider the more general 
map 

x n+ i + x n -i = x n + f(x n ), (10) 

where / does not have to be rational. (Indeed, there is no 
inherent reason to limit the singularity confinement test 
to rational maps, while the notion of degree, and conse- 
quently the definition of the algebraic entropy is tightly 
related to (bi)-rationality. For non-rational transforma- 
tions, one should use a definition of entropy more closely 
inspired from 

It is clear that the map ( |l0| ) can be iterated forwards 
and backwards except possibly when / diverges. Let us 
therefore assume that / diverges at some points x* with 
a power series expansion starting as 

f(x* +e) = a, e -*'(l + 0(e)), Kj > 0, (11) 

and that it vanishes at infinity as 

/(e- 1 ) = 6e L (l + O(e)),L>0. (12) 



From these assumptions it follows in particular that 
f(f(x* + e) + 0(1)) = baT L e LK i{l + 0(e)). The sin- 
gularity analysis proceeds now as follows 



X-l = 


U. 






x = 


Xj + e, 






Xl = 




— u + x*- + e, 




X2 = 


f(x* + < 


- u + baJ L e LK i 


(l + 0(e)), 


X3 = 


- Xj - e - 


^2baJ L e LK >{l + 


0(e)), 


X4 — 


U — Xj + 


A J (e) + 0(e), 





with 

A,(e) = / (-x* - e + 2baj L e LK i(l + 0(e))) - f(x* + e). 

Thus in order for / to pass the singularity confinement 
test, we only need to impose the condition that 

A J (e) = 0(l),Vj, (13) 

and that this term does not cancel the u dependence in 

A simple calculation shows that ( |l3| ) is true at least 
if 1) Kj(L — 1) > 1 and 2) the singularity struc- 
ture of / is even, i.e., both +x* and —x* are singular 
points of / and the expansions at these points match 
as f(Xj + e) - f(-X*j - e) = 0(e). The simplest such 
function is x~ 2 (yielding the map ([j])), but it is easy to 
construct other examples. If attention is restricted to 
rational functions, we can, e.g., pick any two relatively 
prime polynomials Q(x) and P(x) of degree M and de- 
fine f(x) = (P{x)P(-x))/(x 2 Q{x)Q(-x)). As an exam- 
ple we performed the degree growth analysis above on 
the special case P = x + 5, Q = x + 3 and obtained en- 
tropy £ = log ((5 + \/2T)/2). Drawing the orbits again 
corroborates the claim of non-integrability. 

Another class of maps which passes the singularity con- 
finement test is contained in 

X n+ i + X n -i = f n (x n ), (14) 

where we just assume that the functions /„ diverge at 
some points x* (independent of n), and vanishes at in- 
finity, but it is not necessary yet to specify how. Here 
we allow non-autonomous maps, i.e., / may also depend 
on n, as indicated by the subscript (of course this could 
have been done with the previous model as well). The 
singularity analysis now goes as follows: 

X-%= u, 
x =x* + e, 
xi=-u+ f (x* + e), 
x 2 = -Xj - e + fi(fo(x*j + e) - u), 

x 3 =u- f (x* + e) + f 2 (-x* - e + fi(f (x* + e) - u)). 
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For the singularity to be confined at this step, it is suf- 
ficient that the behaviors at the singular points and at 
infinity match so that 

Vj.n lim[/„(z*+e) (15) 

- /„+ 2 (-z* - e + /„+i(/„(x* + e) - u))] = 0, 

The singularity pattern is shorter than for (|l|) , u — > — > 
oo — > — > oo — oo, and upon expanding the ambiguity 
resolves to u. The confinement condition ( |l5|) is similar 
to H). 

The above results show that the singularity confine- 
ment test is only sensitive to the function's behavior at 
its singular points and at infinity. Especially for the non- 
rational case it is easy to dress a function which passes 
the test by something which does not alter this behavior. 

The singularity confinement test is definitely a useful 
tool for identifying potentially integrable systems. It is 
probably necessary but, in the light of the present results, 
it appears to be insufficient. Of course, for a given map, 
the situation would be settled if one could establish any of 
the constructive properties associated with integrability, 
such as Lax pair, superposition principle, conservation 
laws, but in practice this is very difficult. 

It is therefore important to continue developing and 
refining methods for algorithmic testing of integrability. 
For rational maps one such refinement is presented in this 
letter: one should look at the growth of the degree of the 
map (when written in the projective space): if the degree 
grows faster than polynomially (non vanishing "algebraic 
entropy"), it is likely that the system is not integrable. 
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the Institute of Computational Mathematics and Scien- 
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to thank M. Ablowitz, B. Grammaticos, X.B. Hu, S.H. 
Wang and Y.X. Yuan for discussions. 
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